Understanding possible electromagnetic counterparts to loud gravitational wave 
events: Binary black hole effects on electromagnetic fields 
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In addition to producing loud gravitational waves (GW), the dynamics of a binary black hole 
system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plas- 
mas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are 
affected by a pair of orbiting black holes through the merger. In particular, we show how the bi- 
nary's dynamics induce a variability in possible electromagnetically induced emissions as well as 
an enhancement of electromagnetic fields during the late-merge and merger epochs. These time 
dependent features will likely leave their imprint in processes generating detectable emissions and 
can be exploited in the detection of electromagnetic counterparts of gravitational waves. 



I. INTRODUCTION 

The promise of detecting and analysing compact sys- 
tems with both gravitational and electromagnetic waves 
stands out as one of the most exciting prospects in the 
coming decades. As already pointed out in a number of 
works (e.g. [H, 01), most astrophysical systems which pro- 
duce strong gravitational waves likely emit copiously in 
the electromagnetic band. Indeed the strong, and pos- 
sibly highly dynamical, gravitational fields around com- 
pact objects affect the dynamics of plasmas and matter 
which in turn induce different emission mechanisms. 

One example of such a system is a black hole sur- 
rounded by an accretion disk. Strong emission from these 
systems is understood as the result of radiative processes 
within jets powered by the extraction of rotational and 
binding energy. While the latter is qualitatively under- 
stood in terms of Newtonian-rooted arguments based on 
the potential of the central object, the former relies on 
extracting energy from a rotating black hole in the most 
efficient energy convcrtion process we know of. 

The pioneering works of Penrose 0] and Blandford 
and Znajek Q, together with a large body of subse- 
quent work, has provided a basic understanding of pos- 
sible mechanisms to explain highly energetic emissions 
from single black hole systems interacting with surround- 
ing plasmas (e.g. The interaction of electromagnetic 
field lines with the strong gravitational field of a rotating 
black hole is the fundamental component of these mecha- 
nisms to explain the acceleration of particles that traverse 
the black hole's ergosphere. This scenario of a pseudo- 
stationary, single black hole interacting with an accretion 
disk is reasonably well understood, and it is employed to 
explain energetic phenomena such as gamma ray bursts, 
AGNs, quasars, blazars, etc. However, a highly dynam- 
ical stage may occur prior to such a pseudo-stationary 
regime which could give rise to strong emissions. In the 
context of galaxy mergers, such a stage would naturally 



occur as individual black holes in each galaxy eventually 
collide in the galaxy resulting from the merger [||. 

Gravitational waves from such collision would be de- 
tectable by the Laser Interferometric Space Antenna 
(LISA) and, as pointed out in e.g. @, the late or- 
bitting and merger phases would take place within a cir- 
cumbinary disk and possibly even interacting with some 
residual plasma inside the orbiting black holes 

Possible emissions from these systems are only under- 
stood at late times, when the the pseudo-stationary pic- 
ture mentioned above is applicable. However the inter- 
mediate regime has only recently been approached by a 
few works [3, [ToL fll| and in all cases with important sim- 
plifications introduced to track the system. In this work, 
as a follow up of (Io| . we concentrate on understanding 
the dynamics of possible electromagnetic fields anchored 
in the circumbinary disk in the presence of the merging 
black holes. In particular we examine the field configu- 
ration, possible energy enhancement and time variability 
of these fields as the merger take place and point out 
possible process that could give rise to a signal around 
the merger time. To this end, we consider the Einstein- 
Maxwell system in a setup that describes a pair of black 
holes close to the merger epoch, study the electromag- 
netic field behavior and compare to the single black hole 
case. While we do not consider a plasma in our current 
study, our analysis helps to understand the possible be- 
havior in its presence and lay the foundations to future 
work in this direction. 

This work is organized as follows, in Section II, we briefly 
review our formulation and numerical implementation of 
the problem. Section III describes the physical set up to 
study single and binary black hole configurations. Sec- 
tion IV discusses the results obtained for both scenar- 
ios considered and highlight main features which could 
induce emission with particular patterns. We conclude 
with section V which offers some final considerations and 
discussions. 
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II. FORMULATION AND NUMERICAL 
APPROACH 

We solve the coupled Einstein-Maxwell system to 
model the black hole merger interacting with an exter- 
nally sourced magnetic field. We here describe the partic- 
ular formulation of the d equations employed in our sim- 
ulations. We begin with a brief review the 3+1 decom- 
position og Genera Relativity followed by a discussion 
of the Cauchy problem for both Einstein and Maxwell 
equations. 



A. The 3+1 decomposition 

In the Cauchy or 3+1 formulation, the spacetime 
(M, g a b) (a, b = 0, 1, 2, 3) is foliated with spacelike hyper- 
surfaces labeled by constant coordinate time x° = t = 
const. The metric of these hypersurfaces is 7^ = gij 
(i,j = 1,2,3). The normal vector to the hypersurfaces 
is n a = — V a t/| I V a i| I, and coordinates defined on neigh- 
boring hypersurfaces can be related through the lapse 
function, a, and shift vector, /3\ With these definitions, 
the spacetime line element can be expressed as 

ds 2 = g a b dx a dx b 

= -a 2 dt 2 + (dx l + (3 l dt) (dx 3 + ft dt) ; (1) 

while the normal vector/covector is given explicity by 



„« = _ (!,-/?») 
a 



(2) 



Indices on spacetime quantities are raised and lowered 
with the 4-metric, g ab , and its inverse, while the 3-metric 
7ij and its inverse are used to raise and lower indices 
on spatial quantities. The following simple expressions 
relate the 3 + 1 basic variables {"fij, a, f3 1 } with the four- 
dimensional metric {g a b} by 



7 « - 9ij , a = V^Jg™ , ft = Y j g 0j ■ (3) 

In what follows we will make use of both sets of variables. 

B. Einstein equations 

The Einstein equations in the Generalized Harmonic 
formulation [l2|, [l3[ (GH) can be written as a system of 
ten nonlinear partial differential equations for the space- 
time metric g a b- 

g Cd d c d gab + a H b + d b H a = -16 IT \T ab - gab 



+2 T cab H c + 2 g cd g ef d e g ac d f g bd - T ace V 



bdf 



(4) 



where the coordinates x a can be chosen to satisfy the 
generalized harmonic condition 



V c V c x° = -g bc T a bc = H a 



for some arbitrary functions H a . One possibility for de- 
termining these functions, which we adopt here, employs 
the original harmonic condition Hi = (i.e. for the spa- 
tial components), while a damped wave equation for the 
time one Ibi: 



1 



&n c d c H t 



(0) 



Since there is no coupling between the principal part 
of the Einstein equations (|4]) with the generalized har- 
monic condition the full system of equations (Ein- 
stein equations + gauge condition) is trivially hyperbolic, 
determined by wave-like equations with non-linear source 
terms. 

A reduction to first order of the evolution system (ie, 
the Einstein equations (|4|) with the generalized harmonic 
condition © ) can be achieved by introducing new inde- 
pendent variables related to the time and space deriva- 
tives of the fields 

Qab = -n c d c g ab , D iab = d t g a b , (7) 
G = -n c d c H t , Gi = diH t . (8) 

With these definitions we can write the evolution equa- 
tions in our GH formalism in the following way [IH [H[ 



d t g a b 

dtQab 



ft D kab - a Q a b, (9) 

ft d k Qab - af j diD jab 

a d a H b - a d b H a + 2 a T cab H c 

2ag cd D„, !),,>, - Q ca Q db - g ef T ace T bdf ) 



>ld 



Qab - a fW. 



d t D u 



d t H t 
d t G 



- 8ira(2T ab - g ab T) 

- 2a a [n a Zb + n b Z a - g a bn c Z c ] 
+ (7i (3 l (D ia b - d l g ab ), 

= (3 k d k Diab - a dtQab 

+ ^n c n d D icd Q ab + a j jk ?i c D lJC D kab 

- a l a (DM ~ digab), 
= ft G k -aG, 

ft ihC n-.'u),G, 

a f j GiQ jc n c + aT c G c 



(10) 



111) 



{cd 



-Gn c n 
2 

^^(Gi-diHJ-b^^-abG, (12) 

a n 1 



d t Gi 



n a diG 
Gn c n d D lcd + a ^ k G k n c D l]c 



ftd k G 
a 

2 

0-1 a (Gi - diH t ) 



(13) 



This GH formulation includes a number of constraints 
that must be satisfied for consistency. Namely, two sets 
of first order constraints Ci a b,Cijab and the four- vector 
Z a accounting for the physical energy and momentum 
constraints [Ta. [l7|. 



(5) 



c 



dig a b - A 
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Cijab — $i&jab $j&iab 

2Z a ee -T a -H a (t,x i ) = 0. 



(14) 



These constraints are controlled dynamically via the 
inclusion of a constraint damping mechanism [l8j . by 
adding certain terms proportional to these constraints 
(with free parameters <7o and <J\) to the evolution equa- 
tions. 



C. Maxwell equations 

The Maxwell equations can be written in covariant 
form as 



V h F ab = 4TrI a , 



(15) 



\7 b *F ab = (16) 

where F ab is the Maxwell tensor of the electromagnetic 
field, *F ab is the Faraday tensor and I a is the electric cur- 
rent 4- vector. Since F ab is antisymmetric, the divergence 
of equation (|15[) leads to the current conservation 



V Q J Q = 



(17) 



When both the electric and magnetic susceptibility of 
the medium vanish, as in vacuum or in a highly ion- 
ized plasma, the Faraday tensor is simply the dual of the 
Maxwell one, that is 



* j^ab nbed tt< 

r — — 6 r , 



cd 



(18) 



where e abcd is the Levi-Civita pscudotensor of the space- 
time, which can be written in terms of the 4-indiccs Levi- 
Civita symbol rj abcd as 



%bcd 



1 

71 



ibed 



Cabcd = -\fg ff 



bed 



(19) 



For intuitive and practical reasons, it is convenient to 
introduce and work with the electric and magnetic fields 
which are defined as 



E a = F ab n b , B a = *F ab n b 



(20) 



The vectors E a and B a are the electric and magnetic 
fields measured by the normal observer and are purely 
spatial (ic, E a n a = B a n a = 0). The Faraday tensor can 
then be re-expressed as, 

pab = n a E b _ n b E a + fbed ^ ^ ^ ^ 2 1) 
*pab = n a B b _ n b B a _ fbed ^ ^ _ ^ 

We consider here an extended Maxwell system [ijj de- 
fined as 

X7 a (F ab + g ab q>) = -4nl b - a 2 n b ^ , (23) 
= -(T 2 n b (f>, (24) 



V a (*F ab + g ab ' 



which reduces to the standard Maxwell equations if 
"J - = = <p. These extra scalar fields play the role "er- 
ror fields" as they are tightly coupled to constraints vi- 
olations. Moreover, their induced evolution equations, 
obtained by considering the divergence of the extended 
Maxwell equations, imply 



V a V> 



-V„ (<r 2 n a *) 
-V a (<t 2 n a <j>) . 



(25) 
(26) 



which arc a generalization of the telegraph equation. 
Their structure ensures the constraints will propagate at 
the speed of light and will be damped within a timescale 
given by cr^ 1 . This strategy is similar to that defined 
in [H for ideal MHD case. 

The 3+1 version of equations (|23ti24|) . which are the 
ones implemented in our code, are 

(d t -Cp)E i - ,'^v,(n/;, • • ,»-"v,>i< 

atrK E l - 4iraJ i , (27) 
{8 t -Lp)B i + e ijk V j (aE k ) + a 1 ij V j ct> = 

atrKB 1 , (28) 

(d t -Cp)^ + a V i E i = Ana q - aa 2 * , (29) 

(d t -Lp)<t> + a V t B l = -aa 2 <j> . (30) 

with trK the trace of the extrinsic curvature and where 
we have decomposed the current four-vector I a = qn a + 
J a , (with the current J a satisfying J a n a = 0). Obvi- 
ously the standard Maxwell equations in a curved back- 
ground are recovered for '5 = = 0. From equations 
(|2"9")l and (|3"0")) it follows that \I> and <f> can be regarded 
as the normal-time integrals of the standard divergence 
constraints 



ViB 1 = 0. 



(31) 



Finally, just reming that the sources are coupled to the 
geometry by means of the stress energy tensor T ab and 
its trace T = g ab T ab . For our case of interest, the stress 
energy tensor is given by, 



1 

-in 



F n c F b , 



1 



gab F cd F c d 



(32) 



which depends quadratically on the electric and magnetic 
fields. 



D. Implementation 

Our code implements both systems of equations where 
the constraints, as mentioned, are kept under control via 
different but related damping mechanisms: constraint 
damping for the Einstein equations fl8j and the ex- 
tended divergence cleaning for the Maxwell equations as 
explained in the previous section. We adopt boundary 
conditions defined via a combination of Sommerfeld and 
constraint preserving boundary conditions [2l| for both 
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systems. To this end, a characteristic decomposition of 
the Generalized Harmonic formalism eqns. (f9T [T3")) is per- 
formed (at each hypersurface) with respect to the wave 
front propagation direction, given by a normalized spatial 
vector m. This vector is orthogonal to a given boundary 
and belongs to and ordered orthonormal triad {1, p, m}. 
This decomposition gives, 



Dlab , D- 



9ab i Ht 



L a b — Qab — &29ab ± D ma 

L± = G- a 2 H t ± G r , 







(33) 

(i m (34) 
(3 m ± a (35) 
f3 m ± a (36) 



where the symbol {l,p, m} replacing an index means the 
projection along the corresponding vector. The boundary 
conditions are applied only to the incoming modes (i.e., 
{L~ b ,L~} and {D lab , D pab ,G h G p } if (3 m > 0) through 
their time derivatives. We here explain in detail the 
modes related to the metric modes, while the modes re- 
lated to H t are treated in an analogous way. The Som- 
merfeld condition considered is of the type [21| 



(d t + d r + - ){g a b - Vab) = 

r 



(37) 



where rj ab is just the Minkowski metric. By taking a 
time derivative of this equation and rewriting it in terms 
of the incoming characteristic fields, the final form of the 
Sommerfeld conditions results: 



dt 



L ab + (o- 2 - -)g ab 



= 



(38) 



We apply boundary conditions for {Di ab , D pab } if (3 m > 
0. In this case we use the constraint preserving bound- 
ary conditions already given in [Til ], where the original 
time derivatives are corrected by the 4-indcx constraint 
defined in eqn. (|14|) . namely 



(39) 
(40) 



dtDua, = d t D lab - r,„- l>fC, ,, :i , 

d t D pab = dtDp ab - fl m m l Pp[C l3 ab 



where P ab = g ab + n a n b — m a m b is the projection tensor 
on the boundary surface. 

The incoming modes of the electromagnetic fields are 
defined via maximally dissipativc conditions on the time 
derivatives, induced from the physical picture of a cir- 
cumbinary disk present beyond the computational do- 
main. The complete set of eigenvectors for the extended 
Maxwell equations (|27ti30[) . again with respect to ordered 
orthonormal triad {1, p, m}, is given by the following list 
of eigenfields propagating with light speed — (3 m ± a: 



Ei ± Bp 

* ± Err,. 



Ep =F Bi 



(41) 
(42) 



The boundary conditions in this case are just given by 

dtiEi-Bp) = , d t (£ p + B ; ) =0 , (43) 
d t {^-E m ) = , 8t(<f> - B m ) = . (44) 



We adopt Finite Difference techniques on a regular 
Cartesian grid to implement the overall system numer- 
ically. To ensure sufficient resolution is used in an effi- 
cient manner we employ adaptive mesh refinement tech- 
niques. To this end we adopt the HAD computational in- 
frastructure that provides distributed Berger-Oliger style 
Adaptive Mesh Refinement (AMR) [HHl] with full sub- 
cycling in time, together with a novel treatment of arti- 
ficial boundaries [2J|. The refinement regions are deter- 
mined using truncation error estimation via a shadow 
hierarchy [25j and so they adapt dynamically as the evo- 
lution proceeds to guarantee a certain pre-specified tol- 
erance is achieved. A fourth order spatial discretization 
satisfying a summation by parts rule together with a 
third order Rungc-Kutta scheme for the time integra- 
tion are used to help ensure stability of the numerical 
implementation (26[. We adopt a Courant parameter 
of A = 0.2 so that Ai; = 0.2Ax/ at each level= /; 
this ensures the implementation satisfies the Courant- 
Friedrichs-Levy condition dictated by the principal part 
of the equations. However notice that the different damp- 
ing terms (either at the gauge condition ([5]) or constraint 
damping) in turn add a further requirement of the form 
a At i ~ 0(1). Since Ati is considerably larger for the 
coarser grids (where the solution is obtained at large 
distances) we address this issue by making the different 
damping factors space dependent, and in particular 



o l r<r M 

(r-r o M) 2 /(10M) 2 „ ^ „_ Ayr l 40 7 



r > r M 



(i = 0, 1, 2). For similar reasons we also consider a spatial 
dependence for £j. While different options work well, for 
comparison purposes we adopt a similar strategy as in 
thus 



with 



ft = 2 - e is** 



-r 2 /(40M) 2 



(t-t*V 



I _ g (15M) 2 



fi , (46) 



. /a = l; (47) 



with t* chosen to be 10M before the onset of merger. 
As mentioned, other options that "turn on" the gauge 
smoothly from a pure Harmonic condition to a general- 
ized one as given by eqn. $Q work as well. 



III. PHYSICAL SETUP AND INITIAL DATA 

We consider both single and binary black hole simula- 
tions, immersed in an otherwise constant magnetic field 
like the one produced by a disk surrounding the black 
hole at large distances. As the electromagnetic fields in- 
teract with the curved spacetime, they will be dynam- 
ically distorted and eventually reach a quasi-stationary 
configuration. We adopt electromagnetic fields within 
astrophysically relevant values where their energy is sev- 
eral orders of magnitude smaller than the gravitational 
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field energy and so they have a negligible influence on 
the black holes' dynamics. A simple estimate indicates 
this is the case for a large (and certainly astrophysical 
interesting) range of field strengths. To this end first we 
express the field in terms of units of M _1 



Here Bq is the magnitude of the magnetic field at large 
distances from the black hole. For convenience we write 
the solution's components in the standard coordinates 

{t,r,e t <t>h 



B[l/M] = 1.2- 1CT 20 [^L)B[G 

Mq 



(48) 



so, fields up to B Q = 10 18 (M Q /M)G will have energy 
densities < 10~ 4 [7l/~ 2 ]. Additionally, within a sphere of 
radius ~ 100M, the total EM energy will be bounded by 
%1M if the field strength is < 1O 16 (M /M)G. There- 
fore, while we here adopt a field strength of B Q = 
10 4 (M/10 8 Af Q ) G, our results are applicable to much 
stronger values since for fields up to ~ 10 16 (Mq/M)G 
the effects of the electromagnetic fields on the geome- 
try arc negligible however the latter can have a profound 
effect on the former. 

The analysis of the single black hole case will serve 
not only as a test for our numerical implementation, but 
also to understand the features of the initial transient, 
where the EM fields adapt to the geometry of the black 
hole spacetime, giving rise to an electric field and a de- 
formation of the magnetic field (see also [H, 511) ■ As 
it was shown in [28[, the quasi-stationary state is de- 
termined by Wald's solution [3(3 for a Kerr black hole 
immersed in a uniform magnetic field which is aligned 
with its spin. Let us consider the explicit form of Wald's 
solution, which describes a solution of the Maxwell equa- 
tions in the test field case. We assume a Kerr black hole 
in Boycr-Lindquist (BL) coordinates, immersed in a uni- 
form magnetic field [30] . 



ds 2 



2Mr\ 



(r 2 



E 

,2\2 



dt 2 



AM ar sin 2 9 



-dtdcf) 



Aa 2 sin 2 6n 



sin 



A 



dr 2 + me 2 



where E = r 2 + a 2 cos 2 9 and A = r 2 + a 2 — 2Mr. On 
this background, the following tensor defines a solution 
of Maxwell equations: 

F = Fwlj 1 A w° + F 13 uj 1 A lj 3 + F 20 uj 2 A w° 

+ F 23 lo 2 Alu 3 , (50) 

where w°Aw'e \{u a ® uj b - oj h ® w a ), with 



j° = ( — ) (dt - a sin 



Fio 

F13 

F23 



Y}l ,2 d9 



Bo 



ar sm 



sin# 

; EV2 

Ma(r 



[(r 2 + a 2 )d(f> - adt] 



F 



F rt w r Aw' + F r<j) w r AwH F et u e A cj* 



(53) 



F 



in 



A 1 / 2 



F13 , F et = A 1 / 2 F 20 -asin9F 23 , 



F r d> = —a sin 9F 



10 



A 1 / 2 



■ sin 9 F13 
-a 2 )sm9F 



2:5 



Fee = -aA 1 / 2 sin 2 0^0 + (r 2 -( 

Lastly, since we employ excision techniques we prefer 
to adopt horizon penetrating coordinates and so trans- 
form this solution to Kerr-Schild coordinates by the 
transformation 13111 



(54) 



dt 



2r , 

dt —dr , 

A 

dd> -dr . 



(55) 



The explicit expressions of the fields in these coordinates 
are lengthy but straightforward. To gain some insight 
on the solution however, the particular case of a non- 



spinning black hole suffices (setting a = 0). 
element in the BL coordinates reduces to 



The line 



ds 2 



r 2 d9 2 



r J 

2 • 2 
r sm 



dt 2 



r 2 - 2Mr 



dr 2 



(56) 



The electric field vanishes (E l 
(49) field components are given by 



0) and the magnetic 



B r = B a cos 9 , B 6 



-B ^l , ^ 



0. (57) 



The corresponding expressions in Kerr-Schild coordi- 
nates are: 



B r 
E r 



B a cos 9 , B 



E e = , E^ 



■50^,^ = 0,(58) 



2B aM 



(59) 



Notice that although the magnetic field has the 
same expression in both systems of coordinates, there 
appears a toroidal component in the electric field in 
these coordinates due to the non-vanishing shift vector. 



! 0)(1 



cos 



B - 



E 

A 1 ' /2 rsin( 



i*2l 



E 2 

A 1 / 2 asin( 



Bn 



coat 



2 1 2 
r + a 



Bo- 

2Mra 2 (l + cos 



(51) 

To explore the effects of the merger dynamics on the 
electromagnetic field, we compare single spinning black 
hole with cases of equal-mass merging black holes. In 
' all cases, the orbital plane of the evolution (or equato- 
rial plane for the single BH) is assumed to be aligned 
with that of the circumbinary disk. The magnetic field 
is defined as anchored in the disk; hence, its associated 
magnetic dipolc is aligned with the orbital and spin an- 



(52) 



gular momentum. 
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IV. SINGLE BLACK HOLES; THE 
ASYMPTOTIC STATIONARY STATE 

For the single black hole case, we adopt a spinning 
black hole with a spin parameter given by a = 0.7M 
which is close to the spin expected for a merged black 
hole from an equal-mass, non-spinning binary system. 
We adopt this value for comparison with the binary black 
hole scenario presented in the next section, where the fi- 
nal spin can be calculated directly via simulations (for re- 
cent efforts in simulations and data analysis see e.g. [32j 
and references cited therein) or estimated by simple ar- 
guments as in [33|. In this simulation the geometry is 
kept fixed, in order to maintain the same initial choice of 
coordinates and be able to compare easily with Wald's 
(analytical) solution. 

As mentioned, the initial magnetic field is described 
by a poloidal configuration constructed from the electro- 
magnetic potential produced by a circular loop, whose 
radius is assumed to be larger than the region of inter- 
est [Hj]. We assume the disk lies at 10 3 M, and for these 
distances the magnetic field is essentially constant within 
our computational domain, so we simply adopt B l = B a z 
and set the electric field initially to zero throughout . The 
magnetic field strength is B a = 10 G, which is consistent 
with possible values inferred in relevant astrophysical sys- 
tems [35l. |36|. 

Our numerical domain consists of a cubical region de- 
fined by [-80 M, 80 M ] 3 with 80 points in the base grid. 
It employs an FMR configuration with 5 levels of refine- 
ment, each one covering half of the domain of the parent 
coarser level. Thus, the coarsest resolution employed is 
Ax = 2M while the finest one is Ax = 0.I25M. The 
damping parameter is set to be 02 = 1M. 

The evolution shows an initial transient where the 
magnetic fields are deformed through the dynamics ex- 
hibit a twisting behavior around the spinning black hole 
as well as an induced electric field. After t ~ 80M 
the solution is clearly seen to evolve towards a quasi- 
stationary state determined by Wald's solution [30| for 
a Kerr black hole immersed in a uniform magnetic field 
which is aligned with its spin. This is illustrated in Fig.[T] 
which presents both the electric and magnetic field ob- 
tained at t = 200M in the plane y = for x > 0, the 
corresponding field from Wald's exact solution is shown 
for x > 0. The apparent agreement along with a care- 
ful examination of the asymptotic solution indicates that, 
for a black hole immersed within an almost uniform mag- 
netic field aligned with its spin, the final state is Wald's 
solution [H]. 

It is interesting to notice that the electric field configu- 
ration is equivalent to the one produced by a sphere with 
a surface density charge, inmersed in a external mag- 
netic field. This is precisely the result expected within 
the membrane paradigm picture, which endows the event 
horizon with some physical properties (ie, temperature, 
current, density charge, etc.) [371 ] . From this point of 
view, the density charge distribution is thus responsible 




FIG. 1: Magnetic and electric field lines for y — 0, z > at 
t — 200M for a single black hole, together with the apparent 
horizon (green) and the ergosphere (magenta). The regions 
x < display the numerical solution, while Wald's exact so- 
lution is shown in x > 0. 



for the electromagnetic field configuration observed. 



BINARY BLACK HOLES AND EFFECTS ON 
EM FIELDS 



We now turn our attention to the binary black hole 
case. As argued earlier the electromagnetic field con- 
sidered (which is below estimated upper bounds in as- 
trophysical scenarios) is not strong enough to affect the 
dynamics of the black holes, rather the black holes will 
affect the electromagnetic fields. Prior to merger, when 
the black holes are far from each other, the physical pic- 
ture of the EM fields' behavior can be intuitively obtained 
from the knowledge of the black holes quasi-circular tra- 
jectories and the membrane paradigm point of view to- 
gether with results from the previous section. Namely 
the black holes immersed in a uniform field affect the 
electromagnetic field's configuration in its local neigh- 
borhood, and the trajectories of the black holes cause a 
charge separation in the direction perpendicular to both 
the velocity and magnetic field, as in the Hall effect. As 
a result, the electromagnetic field induced by the binary 
motion is given by two dipoles in a quasi-adiabatical 
shrinking orbital behavior. As the merger stage ap- 
proaches, the strong curvature and dynamics might affect 
the EM field's behavior more strongly and so we concen- 
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trate on this stage. We adopt initial data such that the 
merger takes place after about one orbit. This initial 
data corresponds to quasi-equilibrium, equal-mass, non- 
spinning black holes constructed by the publicly available 
LORENE code [38j]. The black holes have masses, given by 
M s = M/2, and are initially separated by « 6M, lying 
beyond the approximate inner most stable circular orbit 
(ISCO) [33} . The initial magnetic field is chosen, as in 
the single black hole case, to be a poloidal configuration 
produced by a circular loop with large radius, so that 
B l = B Q z. The electric field is initially zero through- 
out the computational domain and the magnetic field 
strength adopted is B Q = 10 4 G. 

We adopt a cubical domain given by [—106 M, 106 M} 3 
and employ an AMR configuration with 6 levels of re- 
finement that adjust themselves dynamically to ensure 
that the solution's error is below a pre-determined thresh- 
old using a shadow hierarchy. The coarsest grid has 46 
points, so the coarsest resolution is Ax = 4.6M near the 
boundaries while the finest one is Ax = 0.072M around 
the black holes. We adopt the following set of gauge pa- 
rameters £1 = 0.084/M 2 , £ 2 = 9/M and n = 3, while the 
damping parameters are <7o = o\ = a% = 0.25/M. We 
monitor that the constraint remain well behaved through 
the evolutions; for instance fig. [2] displays the L2 and the 
Loo norms of the physical constraint \\Z\\ (HU), defined 
as 

\\z\\ = Y, z l ■ ( 6 °) 

Notice that besides a small increase at the time of the 
merger and at late times when the finest grid is automat- 
ically discarded and recreated a few times automatically 
as required by the shadow hierarchy, this constraint (as 
well as all others) remains under control during simula- 
tion. 

: 1 1 1 1 1 1 1 1 : 




-80 -40 40 

time(M) 



FIG. 2: Z/2 and the norms of the physical constraint \ \Z\\, 
which is kept under control through the evolution with a small 
increase at the merger. After it, the norms grow and oscillate 
for some period as the finest grid is detroyed and created 
a few times as required by the tolerance adopted at times 
> 38Af the finest is destroyed for the last time and not needed 
anymore. 



A careful inspection of the dynamics described by the 
numerical simulation rcvcalas that, except at late times 
after merger, the physical behavior results significantly 
different than the one corresponding to the single black 
hole case. Indeed the orbiting black holes modify both 
the geometry and the electromagnetic fields. As a result, 
the EM fields are stirred during the evolution, changing 
their configurations as it is displayed in Fig. [3] 




FIG. 3: Magnetic (mostly vertical) and electric field lines 
at different phases during the evolution employing different 
scales for visualization purposes. The figures illustrate dif- 
ferent stages: early when the black holes are separated; near 
merger; shortly after they merge and at late times. The elec- 
tric field lines are twisted around the black hole, while the 
magnetic lines slightly deform from ther initial configuration 
aligned with the z-axis. 

To analyze the influence of the binary's dynamics on 
the electromagnetic field we monitor the (EM) Newman- 
Penrosc radiative quantity $2, and correlate it with ^4, 
which is the gravitational wave Newman-Penrose scalar. 
These scalars are computed by contracting the Maxwell 
and the Weyl tensor respectively, with a suitably defined 
null tetrad 

$2 = F ab n a fh b , * 4 = C abcd n a m b n c m d ; (61) 

extracted at a sphere surface £ located in the wave-zone, 
far away from the sources. We also check that corrections 
required for possible gauge issues, as discussed in (30( | . are 
negligible in our present case. To understand the induced 
multipolar structure of these quantities we decompose 
them in terms of spin- weighted spherical harmonics, with 
spin weight s = — 1 for $2 and s = — 2 for "J 4 (since 
these are their respective spin- weights) . These modes 
exhibit a very similar behavior, with the most relevant 
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FIG. 4: The top figure corresponds to the I — m = 2 modes 
for r^i and r$2 extracted at r — 40M, reescaled properly 
with a factor ~ 3 x 10 6 to fit in the same scale. The bottom 
figure illustrates the phase of these modes. 



ones corresponding to the I = 2, to = ±2 modes, which 
are plotted in fig. [¥] (top). The maximum amplitudes 
of these waveforms correspond to the merger time which 
takes places after over one orbit. Note that since the 
magnetic field is anchored at the disk it does not decay 
with the distance from the binary, which obscures a clean 
interpretation from $2, displaying a non- vanishing m = 
mode at late times, when the stationary state is reached. 
The same happens with the decomposition of the radial 
component of the Poynting vector S r (radial from the 
origin), which shows non-radiative modes not related to 
the binary black hole dynamics. A closer inspection of 
the waveforms (bottom in fig. 2]) reveals that the I = 
to = 2 mode of both the GW and the EM waves oscillate 
with the same frequency, indicating that both are mostly 
dominated by a quadrupolar structure resulting from the 
orbiting behavior. 

The energy carried off by outgoing waves at infinity is 
another interesting quantity. The total energy flux per 
unit solid angle can be found directly from the Newman- 
Pcnrose scalars. 



pGW 

pEM 



dE 



GW 



dtdE 



1 — >oo lt>7T 



^ 4 dt' 



? EM 



(62) 
(63) 



' F™ x factor 

. r- GW 



-40 
time(M) 



FIG. 5: Total flux of gravitational and electromagnetic energy 
corresponding to the binary black hole case (integrated over 
the sphere surface E). The (integrade) electromagnetic energy 
flux has been reescaled by a factor w 10 13 to make it appear 
clearly in the plot. 



located at R% = 30M are shown in fig. [5l both exhibiting 
a maximum at the time of the merger. 

Furthermore, the electromagnetic energy density 
£em = (E 2 + B 2 )/2, while not an invariant quantity in 
general relativity can be employed to get a sense of the 
energy variation in the EM field. As it is illustrated in 
Fig. El as the merger takes place, £em grows significantly 
during the merger and then diminishes as it approaches 
its asymptotic stage which is described by the scenario 
described in the single black hole case. 

Summarizing, several distinct features in the EM fields 
behavior are seen during the evolution: 

• After an initial transient dynamics, and well before 
the merger takes place, the electromagnetic fields display 
a pattern consistent with that produced by equal dipolcs 
orbiting about each other, as it is exhibited in Fig. [7] 
This can be understood from the 'membrane paradigm' 
point of view [37l ] , in which the horizon is endowed with 
a surface density of electrical charge. The quasi-circular 
trajectories of the black holes cause a charge separation 
in the direction perpendicular to both the velocity and 
magnetic field, inducing an electric and magnetic field in 
addition to the external one produced by the far away 
disk. To understand the main features observed a simple 
toy model could be employed. Such model is composed 
of four point charges (A,B,C,D) that orbit circularly. As- 
suming perfect conductivity of the horizon, each charge 
qi(i = A...D) can be written as 



<IA 



-qb = qc 



-Qd, 



(64) 



These quantities, integrated along the sphere surface S, 



where v is orbital velocity of the black holes, B is the 
magnetic field imposed by circumbinary disk and m is 
the radius of the apparent horizon. The orbital motion 
of each charge is written, in Cartesian coordinates, as: 

r'^t') = ifc (cos fii, sin m,0), (65) 

where R A = R Q + r H , R B = Rq — r H , R c = — (i? + 
Th)-, Rd = —{Ro ~ t h), being Rq the distance from 
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FIG. 7: Electric field lines at the plane x = for the binary 
black hole case at different times of the evolution. Notice 
that early on the field configuration is in agreement with the 
expected one from the membrane paradigm with two induced 
dipoles. At the black holes merge the configuration changes 
to a scenario consistent with that required in the Blandford- 
Znajek mechanism. 



FIG. 6: Electromagnetic energy density (normalized to the 
values at the initial time, far distances) at the z — plane for 
the binary black hole case at different times of the evolution. 
Energy contour plots are shown from 1.2 to 0.8 at intervals of 
0.1. 



either black hole to the origin (or center of mass). The 
field produced by this system can be computed by the fol- 
lowing procedure. First, evaluate the Lienard-Wiechert 
potential of each particle and from it obtain its corre- 
sponding electric and magnetic field contribution. Next, 
sum up all contributions from the four particles to ob- 
tain the total electric and magnetic fields. Lastly, add 
the externalt magnetic field to the one produced by the 
charges. The final form of the fields at a given point 
involves implicit equations due to the different contribu- 
tions depending on their respective retarded times. The 
intuitive picture however, can be derived easily by taking 
suitable limits. For instance, when the two black holes 
are well separated (ic, Rq >> rn) the solution around 
each black hole mimics that of a dipolc. At far distances 



from the black holes, both the monopole and dipole con- 
tribution vanish and the system is described to leading 
order by a quadrupole. In this regime, the electric and 
magnetic fields at a point located at ||r|| >> Rq, have 
the dependence 



O 3 

Ew — cos(2Qt + S E ) , B«5„ 
r 



cos(2nt + S B (P6) 



where 8& and 8b are phases whose details are not impor- 
tant for the present argument. This simple form explains 
why fa and 4>4 have the same frequency dependence 2f2. 

• Near and through the merger epoch the fields' 
strengths increase and so does the flux of electromagnetic 
energy. Additionally the fields exhibit a configuration 
consistent with that required by the Blandford-Znajck 
mechanism, namely a (mainly) poloidal/toroidal mag- 
netic/electric fields which could generate interesting en- 
ergy ouput when interacting with a surrounding medium. 
As the merger approaches and take place, the field lines 
become stirred and twisted which, upon rcconncction, 
could also release important amounts of energy, this sce- 
nario is beyond the scope of what we could study with 
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our current setup. 

• Finally, after the merger, the system can be under- 
stood again as a conducting sphere spinning in a external 
magnetic field. Its behavior can be understood by the re- 
sults presented in the previous section. 



VI. CONCLUSION 

We have analyzed the behavior of electromagnetic 
fields influenced by the dynamics of a binary black hole 
system. Our study illustrates several interesting aspects 
of such systems that emit not only gravitational waves, 
but can also radiate electromagnetically. Gravitational 
waves would be emitted through the different dynami- 
cal stages of the system -inspiral, merger and ringdown- 
through a rather smooth manner with the peak occur- 
ring at the merger stage. Electromagnetic waves on the 
other hand, will be emitted through diverse processes 
driven by the interaction of the EM fields with surround- 
ing plasma, gas or matter. In the current work we have 
studied the behavior of the EM fields and illustrated their 
"radiative" behavior as energy that can propagate out- 
wards from the system as the black holes influence them. 
This energy will likely be absorbed and re-emitted by 
the surrounding plasma and would, in turn, be possible 
to observe. The orbiting behavior however, would leave 
its mark in a time variability naturally induced in the 
emission process. 

Indeed, as we have illustrated here, the EM fields have 
a clearly discernible pattern tied to the dynamics of the 
system, making them possible tracers of the spacetime 
-in the electromagnetic sector- as these features would 
imprint particular characteristics in processes producing 
observable EM signals. In particular, in the pre-merger 
stages, the black hole dynamics induce EM flux oscilla- 
tions with a period half that of the dominant GW signal 
produced by the system, i.e., a fourth of the orbital 
period, and a gradual enhancement of the energy in 
the electromagnetic field. This enhancement, together 
with a flux of electromagnetic energy would impact 
surrounding plasma in a stronger way than would be 
the case for a single stationary black hole as the latter 
would neither exhibit an enhancement, nor would it 
give rise to an outward flux of EM energy. Emissions 
in this later case require mechanisms like accretion or 
Blandford-Znajck to take place. Certainly, the same 
requirements will apply to binary black hole systems 
at late times as they give rise (generically) to a single 
spinning black hole. Furthermore, as the merger takes 
place fields are significantly twisted and stirred opening 
up the possibility of interesting emissions through mag- 
netic reconnection. A study of such scenario requires a 
resistive treatment of the problem which is beyond the 
scope of our current work. 

Perhaps even more exciting is the possibility of induc- 
ing a Blandford-Znajck analog for binaries as the merger 



proceeds. As we have seen, the system's dynamics 
induces a configuration consistent with the basic picture 
of this process and, by extracting rotational energy 
from the system, more powerful emissions could be 
expected. For this to take place however, an important 
conservative requirement should hold. Namely, an 
crgosphcre must be present so that rotational energy 
can be extracted from the merging black holes. The best 
case scenario is for the ergosphere to form before the 
plunging phase begins so that the black holes can orbit 
and sufficient time is left for the extraction to occur. An 
estimate for when this takes places can be drawn from 
an "effective one body" approach, as described in the 
appendix, which indicates highly spinning configurations 
are required, at least in the particle limit, for an orbiting 
behavior to exist within the ergosphere. 

Finally, and at a rather academically interesting level, 
it is interesting to ask what conditions would be required 
to extract so much energy that the end state is essentially 
a non-spinning black hole after the merger. For this to 
take place, the timescale of the BZ process (tbz) should 
be comparable to that of the merger (tm)- The latter 
typically lasts tm — 50M, the former can be estimated 
as in [f| giving rise to 

T BZ 5 x lO 8 (lO 15 G/B) 2 (Af /M)M , 
~ lO 7 (lO 15 G/B) 2 (M /M) 2 r M ; 

thus, for B > 10 10 G both times result comparable for BH 
masses > lO 8 Af . Alternatively, one can estimate the 
amount of energy that could be extracted as a function 
of field strength within tm = 50M. If M p is the final 
irreducible mass of the black hole after rotational energy 
has been extracted via the BZ mechanism, M p obeys 



M, 



M - M, 



10 8 (10 15 G/B) 2 (M Q /M) 5 



(67) 



Thus, for M = 1O 9 M and B = {10 7 /10 8 }G about 
{10- 4 /10- 2 }M o c 2 (i.e. ~ {10 48 /10 50 } ergs) is released 
in roughly a day. 

While fields of these magnitudes might be unlikely, it 
is interesting that the strengths are not completely out 
of nature's ability to manifest. Beyond this possibility 
the impact of the dynamics on the electromagnetic fields 
is to induce a distinct variability as fields are dragged 
by the black holes. This suggests tantalizing prospects 
to detect pre-merger electromagnetic signals from sys- 
tems detectable in the gravitational wave band. How- 
ever, a complete description of the problem requires the 
incorporation of gas and radiation effects. Notwithstand- 
ing these missing -and important- ingredients, the main 
qualitative features — driven by the orbiting behavior of 
the black holes, whose inertia is many orders of magni- 
tude above all else would intuitively remain unaltered. 

At a more speculative level, these combined signals 
could be exploited to shed light on possible observations 
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to analyse alternative theories of gravity where photons 
and gravitons might propagate at different speeds or 
gravitational energy could propagate out of our possi- 
ble 4-dimcnsional brane (for a recent discussion of some 
possibilities see |40ll4lj|). 

As a final comment we stress that while our work 
is a step towards understanding possible emissions 
induced by binary black hole merger processes, we have 
only scratched the surface of possible phenomenology. 
For instance, within the current approach, scenarios 
with unequal masses and/or spins must be investigated 
and work is in progress to address them [42j. Still, a 
complete understanding of associated phenomena will 
require investigating, in particular, the interaction with 
a surrounding plasma and associated possible emission 
mechanisms. A preliminary related step in this direction 
has been considered in [{J [43[ where the possibility of 
emission by a fossil gas in between the black holes has 
been considered. 
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VII. APPENDIX 

The location for the innermost stable circular orbit, at 
the equator, in a black hole spacetime of mass M and 
spin parameter a is given by [44J 

z ^ l+ { l -W) [( 1 + m) + ( 1 ~m) 




risco = il/(3 + Z 2 T[(3-Zi)(3 + Zi+2Z 2 )] 1 / 2 ) 
The ergosphere, on the other hand, is located at 



tergo =M+ ^AP - a 2 cos(#) 2 

Thus, at the equatorial plane tergo = 2M while for 
prograde orbits rjsco( a — 0) = 6A/ and rjsco(a- — 
M) = M. Consequently, for sufficiently high spins 
tergo > risco- This is illustrated in figure [8j where 
the critical value at which the two lines cross A related 
interesting point is that this argument bears relevance 
also to the generation of gravitational waves themselves. 
If the black holes orbit outside the isco but inside the er- 
gosphere they could tap rotational energy and produced 
stronger emissions. All cases so far studied numerically 
lie below this critical value and so, if this simplistic model 
holds, it would indicate simulations have not yet probed 
the possibility of extraction of rotational energy. There- 
fore binary black hole systems could still potentially yield 
further interesting features in such regime. 




a/M 

FIG. 8: ISCO and ergosphere radii vs a/M. For a > 0.943M 
the ergosphere lies beyond the ISCO. 
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